arXiv:1502.03992vl [cond-mat.mtrl-sci] 13 Feb 2015 


Stick-Slip Control in Nanoscale Boundary Lubrication by Surface Wettability 
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We study the effect of atomic scale surface-lubricant interactions on nanoscale boundary- 
lubricated friction, by considering two example surfaces - hydrophilic mica and hydrophobic 
graphene - confining thin layers of water in molecular dynamics simulations. We observe stick- 
slip dynamics for thin water dims confined by mica sheets, involving periodic breaking-reforming 
transitions of atomic scale capillary water bridges formed around the potassium ions of mica. How¬ 
ever, only smooth sliding without stick-slip events is observed for water confined by graphene, as 
well as for thicker water layers confined by mica. Thus, our results illustrate how atomic scale details 
affect the wettability of the confining surfaces, and consequently control the presence or absence of 
stick-slip dynamics in nanoscale friction. 

PACS numbers: 62.20.Qp, 68.35.Af, 68.08.Be 


Understanding friction plays a central role in tech¬ 
nological applications and phenomena in diverse fields 
ranging from micromechanical devices to bioengineering 
[l[ and to earthquakes [2|. Given the continuing minia¬ 
turization of mechanical devices towards the nanoscale 
3, improved understanding of friction and wear could 
help in reducing energy consumption, improving reliabil¬ 
ity and extending service life. Indeed, an important part 
of their design process consists of trying to minimize fric¬ 
tion and to eliminate stick-slip dynamics [^. 


Stick-slip control in lubricated friction is of particular 
importance given the vast amount of applications where 
lubricants are used to reduce the detrimental effects of 
friction and wear [^. Examples of mechanisms behind 
the emergence of stick-slip in boundary-lubricated sys¬ 
tems have been numerically demonstrated to include re¬ 
peated crystallization and shear melting of the thin lu¬ 
bricant film 0, interlayer slips within the ordered solid¬ 
like lubricant film, or wall slips at the wall-film interface 
0] . Most of the numerical studies of stick-slip in bound¬ 
ary lubrication have focused on coarse-grained or simpli¬ 
fied/idealized models [1, explicitly considering 

the atomic scale interactions occurring in real systems. 
On a coarse-grained scale, a useful classification of the 
lubricant-surface interactions is given by the wettability 
of the confining surfaces by the lubricant, with systems 
displaying a larger contact angle/lower wetting gener¬ 
ally exhibiting lower friction. Other approaches to fric¬ 
tion control include e.g. applying mechanical oscillations 
[IS El- While the effect of wettability on lubricated 
friction has been studied experimentally in macroscopic 
[iS - flsjl and nanoscale [lS| systems, and modeled using 
phenomenological finite-element models 17| and simpli¬ 
fied molecular dynamics (MD) simulations of nanopat- 
terned surfaces [isl ll9| , less is known about the under¬ 
lying atomic scale processes and mechanisms responsible 
for the presence or absence of stick-slip. 


Given the large surface-to-volume ratio in boundary lu¬ 
brication, nature of the interaction between the lubricant 
and the confining surfaces originating from their atomic 
composition should play a crucial role. Thus, we study 
the interaction of a thin water layer (thickness h around 
0.5 nm unless stated otherwise) in MD simulations us¬ 
ing full atomic models of two experimentally relevant 
confining surfaces with different wetting characteristics: 
crystalline miem a hydrophilic substrate that strongly 
adsorbs water [20| and graphene, a hydrophobic surface 
interacting weakly with water 2l|, see Fig. [TJ We ob- 



FIG. 1. (color online) (a) The geometry of the simulation 
system. Solid sheets are held together by a constant normal 
load Fn. The top sheet is moving at a constant velocity V, 
and the bottom sheet is connected to a fixed stage by a spring 
of stiffness k. The water molecules are confined by (b) two 
mica sheets (each of thickness of 8.34 A) or (c) two monolayer 
graphene sheets. The color code of the atoms is: water oxygen 
(red), water hydrogen (white), potassium (pink), silicon (yel¬ 
low), aluminum (blue), mica oxygen (cyan), mica hydrogen 
(lime), and carbon (gray). 
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serve stick-slip dynamics for thin water layers confined 
by mica: each unit cell of mica contains two K"*" ions, 
interacting strongly with the water oxygens via Coulomb 
interactions, leading to formation of atomic scale capil¬ 
lary bridges next to the K’*' ions, connecting the two mica 
surfaces in the stick state. These bridges break during 
the subsequent slip event, and reform during the next 
stick phase, a process that is also visible as the breaking 
and reforming of interfacial hydrogen bonds between wa¬ 
ter and mica. This mechanism is different from both the 
crystallization-shear melting transitions @ and interlayer 
or lubricant-surface slips Q observed before in simplified 
models. In contrast, water films confined by hydropho¬ 
bic graphene, as well as thicker water layers confined by 
mica, exhibit fundamentally different dynamics with no 
stick-slip. 


To model the conhned water him, we consider sys¬ 
tems ranging from 200 to 1200 SPC/Fw water molecules 
(^ . We consider 2Ml-muscovite mica with the formula 
KAl 2 (Al, S^Oio(OH) 2 , with the force held parameters 


from Ref. 


One mica surface consists of 10 x 6 unit 


cells, and has linear dimensions of Lx = 52.07 A and 
Ly = 54.036 A, see Fig. [TJ To create site disorder, 
mimicking a real mica surface with a random distribu¬ 
tion of potassium ions on it, one K+ ion of the pair in 
each unit cell is removed and subsequently placed on the 
bottom part of the sheet [1^. The graphene sheets have 
Lx = 68.063 A and Ly = 36.841 A. The Lennard-Jones 
parameters for carbon are from Ref. 2^. The cutoff ra¬ 
dius is Tc = 10.0 A for all potentials. Both sheets are par¬ 
allel to the xy plane with periodic boundary conditions 
along the x and y directions. Couette how is generated by 
moving the top sheet at a constant velocity V along the x 
direction. The distance between parallel sheets is allowed 
to vary, and a constant normal load Fk, giving rise to a 
pressure P±, is applied on the top sheet. The bottom 
sheet is constrained to move along the x axis, and is at¬ 
tached to a spring of stiffness k/Np = 0.0035 N/m, where 
Np is the total number of atoms in a sheet. The other end 
of the spring is connected to a hxed stage. Temperature 
of T = 295 K is maintained using a Langevin thermo¬ 
stat, app^d only in the y direction to avoid streaming 


bias |26l.l27|. The equations of motion are solved with the 


velocity Verlet algorithm implemented in the LAMMPS 
code 2^, with an integration time step of 1 fs. Long- 
range electrostatic interactions are computed using the 
particle-particle particle-mesh (PPPM) solver with 10“® 
accuracy. Initially the water molecules are arranged in 
a simple cubic lattice. The simulations are first run for 
100 ps with both surfaces kept fixed, followed by 100 ps 
during which the top surface is subject to a normal force 
An and is allowed to move vertically. Then, the top sur¬ 
face is driven horizontally with a velocity V for 1 ns to 
generate the steady state, after which we continue the 
simulations for approximately 60 ns, recording the ob¬ 
servables of interest. 
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FIG. 2. (color online) Time evolution of (a) the spring force 
per sheet atom, (b) the friction force per sheet atom on the 
bottom mica sheet applied by the water and the top mica 
sheet, (c) the position Z of the center of mass of the top 
mica sheet in the ^ direction, and (d) the number of hydrogen 
bonds between the 256 mica-confined water molecules and the 
bottom mica sheet, (e) Spring force per sheet atom for 200 
water molecules confined by graphene sheets. F = 0.1 m/s 
and P± = 1 atm in both cases. 


Simulation results for 256 water molecules confined by 
mica sheets for P± = 1 atm and F = 0.1 m/s are shown 
in Fig. O The force per atom on the bottom sheet ap¬ 
plied by the spring, F^/Np, exhibits characteristic stick- 
slip behavior [Fig. [2] (a)]. Fig. [2] (b) shows the friction 
force per sheet atom on the bottom mica plate applied by 
the water and the top mica plate, Fj-fNp, exhibiting sim¬ 
ilar time-dependence as the spring force, with superim¬ 
posed high-frequency fluctuations due to the finite tem¬ 
perature. Fig. [2](c) shows the position Z of the center of 
mass of the top sheet in the z direction. The center of the 
bottom mica sheet is fixed at z = 4.16 A. During each 
slip event, Z increases by roughly 10% Q. Since forma¬ 
tion and breaking of interfacial chemical bonds is known 
to play a role in friction (see Ref. [2^ for an example 
from rock friction), we show also the time-dependence of 
the number of hydrogen bonds (i.e. the number of water 
hydrogens closer than 3 A from the bottom mica surface) 
between water and the bottom mica surface in Fig. [11(d): 
bonds break as the system evolves from stick to the slip 
state. 

For comparison, we also performed MD simulations of 
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water confined by hydrophobic graphene sheets. We var¬ 
ied the number of water molecules from 200 to 1200, the 
normal loads from = 1 to 10 atm, and the driving 
velocities from V = 0.01 to 0.1 m/s. Fig. [5] (e) shows 
the spring force from simulations of 200 water molecules, 
P± = 1 atm, and 1^ = 0.1 m/s; similar results are ob¬ 
tained for other P± and V values. We observe a small 
increase of friction with V for both mica and graphene, 
see Supplemental Material [ ' 


and Refs. 31, M for 


experimental results on mica-confined systems with slid¬ 
ing velocities significantly lower than those reachable in 
our MD simulations. Our simulations thus demonstrate 
that the stick-slip behavior does not arise for thin water 
films conhned by graphene. Instead, continuous, smooth 
sliding with the maximum friction force well below that 
obtained for mica is observed for all parameter values 
considered. We also note that the same applies to the 
mixed system with one graphene and one mica surface: 
slip is localized at the hydrophobic graphene-water inter¬ 
face, and no stick-slip is observed. 

This difference between the two kinds of surfaces may 
be explained by the relatively strong interaction of the 
potassium ions on the mica surfaces with the oxygen 
atoms of the water molecules via Coulomb interactions. 
Thus, the ions could act as “freezing nuclei”, with the wa¬ 
ter molecules gathering around them to form nanoscale 
capillary water bridges [H, , connecting the top and 

bottom surfaces within the stick phase. As the system 
starts to slip, these bridges would break. The interac¬ 
tion of carbon atoms with oxygen is much weaker, and 
we expect that no capillary bridges are formed between 
graphene sheets, explaining the absence of stick-slip dy¬ 
namics in that case. 

To verify this hypothesis, we calculate the density dis¬ 
tributions p(x, y) of water oxygens in the contact layer 
relative to the bottom surfaces. Fig. [3] (a) shows p{x, y) 
for a water film confined by mica sheets when the system 
sticks [f = 1 ns in Fig. [2] (a)]. Peaks in p{x, y) are located 
at the K+ ions. Fig. [3] (b) presents the corresponding 
p{x, y) during the first slip state when t = 5 ns [cf. again 
Fig. [2] (a)]: the peaks of p{x,y) become smaller and 
broader. Finally, Fig. (c) shows p(x, y) for the subse¬ 
quent stick state at t = 7 ns [Fig. [2] (a)], where we again 
observe that the peaks are as high and narrow as those 
of the previous stick state. 

To gain more insight into the nucleation and breaking 
of the capillary bridges between the surfaces, we calcu¬ 
late the density profiles p{z) of water oxygens across the 
gap. When the system is slipping [Fig. |3](d)], p{z) ex¬ 
hibits two separate peaks, consistent with breaking of 
the capillary bridges. In the stick state [Fig. [3] (e)], 
p{z) exhibits multiple peaks spanning the gap. This can 
be understood as the water molecules forming nanoscale 
capillary bridges between the two mica surfaces. In con¬ 
trast to this behavior, the density distributions p{x, y) of 
the water film consisting of 200 water molecules confined 





FIG. 3. (color online) Contour graphs of the density distri¬ 
bution p{x, y) of water oxygens in the contact layer relative 
to the bottom mica surface for (a) t = 1 ns (“stick”), (b) 
t = 5 ns (“slip”), and (c) t = 7 ns (“stick”). White cor¬ 
responds to no water molecules being present. The density 
profiles across the gap, p{z), of water confined by mica sheets 
when the system (d) slips and (e) is in the stick state. In 
both (d) and (e), the top surface of the bottom mica sheet is 
a.t z = 8.3 A, while the lower surface of the top mica sheet is 
a.t z = 14.7 A in (d) and a.t z = 13.4 A in (e). 


by graphene sheets in Fig. 0] show that water clusters 
to form a single, relatively large droplet-like structure 
between the two graphene sheets, without any apparent 
signature of breaking-reformin g tr ansitions. The corre¬ 
sponding density profiles p{z) (30| are similar to previ¬ 
ous observations in equilibrium graphene-confined sys¬ 
tems [ 3 ^. 


Thus, when the two mica surfaces are very close to¬ 
gether, the thin confined water film loses its fluidity, and 
the bulk flow properties of water play little or no role 
in friction. However, they may be recovered by increas¬ 
ing the thickness of the water layer 36|, with the con¬ 


ditions approaching those of hydrodynamic lubrication. 
To this end, we performed MD simulations with four dif¬ 
ferent, larger thicknesses of the water layer: h = 1.77, 
2.03, 2.29, and 2.56 nm, corresponding to 1536, 1792, 
2048, and 2304 water molecules, respectively. For these 
thicker water films, the stick-slip dynamics disappears. 
Instead, smooth sliding dynamics is observed, which at a 
first glance looks similar to that in the graphene-confined 
system. However, subtle differences can still be observed 
between the two surfaces. Zooming in to the spring 
force time series (e.g. the one shown in Fig. IKe)] re¬ 
veals periodic oscillations corresponding to the eigenfre- 
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FIG. 4. (color online) Contour graphs of the density distri¬ 
bution p{x, y) of water oxygens in the contact layer relative 
to the bottom graphene surface for (a) t = 0 ns, (b) t — 3 ns, 
(c) t = 5 ns and (c) t = 8 ns. White corresponds to no water 
molecules being present. 


quencies of the spring-bottom plate mass (M) system, 
/ = l/(27r) For both surfaces, the varying am¬ 

plitudes of these oscillations at each period [blue circles 
in Figs. [5l[a) and (b)] form sequences of time-ordered 
observations X (n) which can be well-described by an au¬ 
toregressive model X{n + 1) = aX (n) + W (n) [or equiva¬ 
lently, the Ornstein-Uhlenbeck process X{n+1)—X{n) = 
— {1 — a)Xin) + W], with W white noise originating from 
the interaction with the fluctuating lubricant and a a 
model parameter, both extracted using the R package 
[s^- For water confined by graphene, we find a ~ 0.8 
and 5W « 0.1 pN for all conditions considered, while 
we find a k, 0.1 and 6W « 0.3 pN for thick water films 
[h > 1.77 nm) confined by mica. Accordingly, the auto¬ 
correlation function (ACF) of X{n) for mica decays more 
rapidly to zero than its counterpart for graphene. In both 
cases the ACFs computed from the simulation data agree 
with those of the corresponding autoregressive model [see 
Figs. [5](c) and (d)]. The observation that 5W does not 
significantly depend on h for h > 1.77 nm indicates that 
the screened mica-water interaction has a sub-nanometer 
range, resulting essentially in a surface effect of the fluc¬ 
tuations of the water layer. Also, the stronger interaction 
of mica with the fluctuating lubricant (as compared to 
that of graphene) results in a factor of three larger 6W 
[see also Figs. [5](e) and (f)]. 

In summary, the presence or absence of breaking¬ 
reforming transitions of local capillary bridges in the 
water film, controlled by the atomic structure and the 
ensuing wettability (hydrophilic mica vs hydrophobic 
graphene) of the confining surfaces, plays a crucial role 
in whether stick-slip dynamics is observed or not. For 
mica, the decisive role of the K+ ions in the formation of 
the nanoscale capillary bridges suggest that the micro¬ 
scopic details behind stick-slip dynamics should in gen¬ 
eral depend on the atomic structure of the system, and it 



FIG. 5. (color online) Time dependence of Fa/Np during a 
time period of 0.1 ns (a) for 1200 graphene-confined and (b) 
for 1536 mica-confined water molecules. Blue circles show 
the local maxima of the signals, corresponding to the time- 
varying amplitudes X (n) of the spring force oscillations. The 
autocorrelation functions (as a function of the lag t = n — n') 
of these amplitudes are given (c) for graphene and (d) for 
mica, extracted from 10 ns long spring force signals. The solid 
lines correspond to the simulation results, while the dashed 
lines show the corresponding AGFs from the autoregressive 
model. The plots of X[n + 1] vs X[n\ extracted from the 
simulations for (e) graphene and (f) mica further illustrate 
the different nature of the smooth sliding dynamics for the 
two kinds of confining surfaces. The slopes of the lines (linear 
fits) are 0.9 and 0.1 for graphene and mica, respectively. 


would be interesting to perform similar studies for other 
confining surfaces with different surface-lubricant inter¬ 
actions. Nevertheless, we expect our irrain observations 
to be rather general, and to open up interesting possibili¬ 
ties in controlling nanoscale boundary-lubricated friction 
by tuning the wettability of the confining surfaces. 
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